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We analyze the two-dimensional momentum distribution of electrons ionized by few-cycle laser 
pulses in the transition regime from multiphoton absorption to tunneling by solving the time- 
dependent Schrodinger equation and by a classical-trajectory Monte Carlo simulation with tunneling 
(CTMC-T). We find a complex two-dimensional interference pattern that resembles ATI rings at 
higher energies and displays Ramsauer-Townsend diffraction oscillations in the angular distribu- 
tion near threshold. CTMC-T calculations provide a semiclassical explanation for the dominance 
of selected partial waves. While the present calculation pertains to hydrogen, we find surprising 
qualitative agreement with recent experimental data for rare gases Q. 

PACS numbers: 32.80.Rm,32.80.Fb,03.65.Sq 



The interaction of few-cycle laser pulses with matter 
has recently attracted considerable interest Q as increas- 
ingly shorter pulses with duration of the order of 10 fs and 
below became available. Novel aspects of laser-matter 
interactions such as the dependence of high-harmonic ra- 
diation or electron emission on the carrier envelope phase 
0, 0] and the interference of electronic wave packets 
emitted at different points in time during the ultrashort 
pulse became apparent. Another recent advance is 
the imaging of the momentum distribution of the ionized 
electron providing insight into the ejection of both one- 
electron [lj and non-sequential multiple electron emis- 
sion 0. For single-electron emission, the longitudinal 
momentum distribution (fc z along the direction of the 
laser polarization) of photoelectrons from rare gases fea- 
tures a broad "double-hump" structure near threshold 
which surprisingly resembles the k z distribution for non- 
sequential double ionization Q • While for the latter case 
this structure results from electron-electron collision dur- 
ing rescattering of the laser driven electron at the ionic 
core, in the former case it is due to the interplay of the 
Coulomb interaction and laser field on the receding tra- 
jectory 0,H,SE3- By contrast, the transverse momen- 
tum distribution features a narrow Coulomb-like cusp 
well-known from ion-atom collisions 
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recently, Rudenko et al. [l| presented first fully two- 
dimensional (k z ,k p ) momentum maps for laser- ionized 
electrons from different rare gases, displaying a complex 
pattern whose origin is, so far, unexplained. In this letter 
we investigate the 2D momentum map for laser ionization 
of hydrogen and find an equally complex yet surprisingly 
similar pattern suggesting a simple explanation indepen- 
dent from specific properties of the target atom. 

The interaction between the laser field and the atom 
can be characterized by two different mechanisms con- 
trolled by the value of the Keldysh parameter 7 = 
\J I p /2Up, where I p is the ionization potential of the 
atom, Up — Fq/Auj 2 the ponderomotive energy, oj the 



laser angular frequency, and Fq the peak amplitude of 
the laser field. In the multi-photon regime (7 > 1) the 
atomic interaction is governed by the quantum nature of 
the radiation field resulting in the absorption of n pho- 
tons (n > 1) from the field. By contrast, in the tunneling 
regime (7 < 1) the atom responds to the strong per- 
turbation by a "classical" electric field where ionization 
proceeds via tunneling. The present calculation as well 
as recent experiments explores the transition regime 
around 7 w 1 where a more complex response is to be 
expected. We focus on hydrogen in order to avoid any 
ambiguity resulting from additional approximations re- 
quired for many-electron targets. 

The Hamiltonian of a hydrogen atom driven by a lin- 
early polarized laser field is 

H = \ + V{r) + zF{t), (1) 

where V(r) = —1/r is the Coulomb potential energy, 
p and f are the momentum and position of the electron, 
respectively, and F(t) is the time dependent external field 
linearly polarized along the z direction. The laser pulse 
is chosen to be of the form 

F(t) = F sin 2 ^— J cos(ut + tp) (0 < t < r) , (2) 

where u> is the laser frequency, ip = the carrier-envelope 
phase, r the total pulse duration, and Fq is the peak field. 
Atomic units are used throughout. 

The time-dependent Schrodinger equation (TDSE) can 
be solved by different techniques |l.T , ll4j . Approximation 
methods include semiclassical approximation methods 
[H E3, the (Coulomb-) Volkov approximation [HI E3] 
and CTMC-T method HE^. We employ the generalized 
pseudo-spectral method for solving the TDSE [20] . The 
method combines a discretization of the radial coordi- 
nate optimized for the Coulomb singularity with quadra- 
ture methods to allow stable long-time evolution using 
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a split-operator method. Both the unbound as well as 
the bound parts of the wave function \tp(t)) can be accu- 
rately represented. The calculation of the 2D momentum 
distribution requires projection of the partial waves |fc, 1} 
onto outgoing Coulomb waves 
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after the conclusion of the pulse. In Eq. Si(k) is 
the momentum-dependent Coulomb phase shift, 9 k is the 
angle between k and the polarization direction of the 
laser field, z, Pi is the Legendre polynomial of degree 
I, and \k,l) is the eigenstate of the atomic Hamiltonian 
with positive eigenenergy E = fc 2 /2 and orbital quan- 
tum number I. The atom is initially in its ground state. 
Due to the cylindrical symmetry of our system for a lin- 
early polarized laser field, m — is a constant of motion 
and omitted. The distortion of the momentum distribu- 
tion due to long-range final-state Coulomb interactions 
is fully accounted for Eq. ©. We also have performed 
classical-trajectory Monte-Carlo calculations 9j incorpo- 
rating tunneling (CTMC-T) which include both Coulomb 
and laser field interaction non-perturbatively. 

Examples of the two-dimensional momentum distribu- 
tion {k p ,k z ), d fj kz = 2nk p (^P) , for an 8-cycle pulse 

(t = 1005), frequency Li = 0.05 and different field am- 
plitudes F = 0.0377(7 = 1.34), 0.0533 (7 = 0.95), and 
0.075(7 = 0.67) are shown in Fig. 1, illustrating the 
transition from the multi-photon to the tunneling regime. 
Each frame displays a complex interference pattern which 
is characterized by a transition from a ring-shaped pat- 
tern at larger k = \J k 2 + k 2 > 0.4 with circular nodal 

lines to a very different pattern of pronounced radial 
nodal lines for small k near threshold. The first point 
to be noted is that the overall pattern displays a surpris- 
ing and striking similarity to the experimental pattern 
observed recently for rare gases, such as helium, neon, 
and argon The ring pattern is reminiscent of ATI 
peaks of the multi-photon regime. The point to be no- 
ticed is that ATI rings are present well into the tunnel- 
ing regime for larger energies E = ifc 2 > 1. Concurrent 
CTMC-T calculations for the present laser parameters 
show that classical rescattering provides only a minor 
contribution to the electron spectra above E = 2U p . En- 
ergies E >2U P are effectively inaccessible by the classical 
quiver motion thus rendering the multiphoton quantum 
process as dominant pathway. The ATI-like component 
can be quantified by the partial ionization probability p\ 
with orbital quantum number I residing within a given 
i-th ring between adjacent minima (ki — A.;, ki + Aj) 
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The i-th circle has a mean radius fcj = y/2Ei, where the 
energy Ei corresponds to the i-th ATI peak of the pho- 
toelectron spectrum. For i > 2 (k > 0.4) rings become 
recognizable. 
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FIG. 1: Doubly-differential electron momentum dis- 
tributions in cylindrical coordinates (k z ,k p ). The 
parameters of the field are u> — 0.05, t = 1005. In (a) 
7 = 1.34 (F = 0.0377), (b) 7 = 0.95 (F = 0.053), and 
(c) 7 = 0.67 (F = 0.075). 

The transition to an entirely different and unexpected 
radial pattern occurs for energies within « 2lo from the 
threshold. The radial nodal pattern at low energies can 
be made more explicit by analyzing the angular differ- 
ential probability, d 2 P/dkd(cos9k), at fixed k (Fig. 2). 
The probability displays pronounced oscillations that re- 
markably resemble those of a single Legendre polynomial, 
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At lower field and higher k = 0.34, the Legendre polyno- 
mial with Iq = 6 dominates while at higher field and 
low k = 0.19 the polynomial lo — 8 has the largest 
weight. The dominance of a single Legendre polyno- 
mial Pi implies the dominance of a single or a few par- 
tial waves in the momentum-differential ionization cross 
section. The resulting radial nodal pattern with pro- 
nounced minima at certain angles is a well-known feature 
in low-to-intermediate energy electron-atom scattering 
referred to as generalized Ramsauer-Townsend (CRT) 
diffraction oscillations 0, H3, Hjl . The present result 
suggests that laser-driven scattering of electrons in the 
Coulomb field of the nucleus leads to similar Ramsauer- 
Townsend-like interference fringes in the angular distri- 
bution. The distribution of contributing partial waves, 
pi, is presented in Fig. 3. Marginally into the multi- 
photon regime (7 = 1.34, Fig. 3a), the I distribution of 
the first ring displays a broad distribution peaking at odd 
angular momenta I = 3, 5, and 7 while in the second ring 
the distribution features a single peak at 1$ = 6 in agree- 
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ment with d 2 P/dkd(cos9k) (Fig. 2a). In the tunneling 
regime (7 = 0.67, Fig. 3b), the first ring near threshold 
peaks at lo = 8. The dominance of a single I is further 
enhanced by the relative suppression of adjacent angular 
momentum Iq±1 of opposite parity, which is a remnant of 
the multi-photon parity selection rule. The second ring 
shows a peak at lo = 9, even though no dominance of a 
single partial wave is evident when looking at the angular 
distribution (not shown). 
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FIG. 2: Angular distribution calculated by solving the 
TDSE (thick solid lines) for u = 0.05, r = 1005. In (a) 
7 = 1.34 (Fq = 0.037) and k = 0.34, and in (b) 7 = 0.67 
(Fq = 0.075) and k = 0.19. Dashed lines: square of 
the Legendre polynomial Pi (cos9k)Jo = 6 in (a), and 
h = 8 in (b). 

GRT interference fringes in electron-atom scattering 
can be semiclassically described in terms of interferences 
of paths with (in general) different angular momenta 
scattered into the same angle 9k fUJ ■ In the special case 
that a single partial wave and thus a single Legendre 
polynomial Pi (cos9k) dominates, the angular distribu- 
tion of the interfering paths must have very similar clas- 
sical "impact parameters" such that they belong to the 
same angular momentum quantum bin [lo, io + l] UM- An 
analogous path interference occurs in laser-atom ioniza- 
tion in the tunneling regime near threshold. To uncover 
the relevant classical paths we employ a CTMC-T sim- 
ulation for the same parameters as in Fig. 3b. The 
ensemble of ionized electrons near threshold features, in- 
deed, an I distribution (Fig. 4a) that resembles the quan- 
tum distribution (Fig. 3b) with a peak near Iq = 8; 
clearly emphasizing the underlying classical character of 
this process. A typical electron trajectory after tunneling 



shows a quiver motion along the polarization of the laser 
field. An interesting observation is that, even though the 
motion is strongly driven by the laser field, the motion 
follows the Kepler hyperbola (Fig. 4b). The point to 
be emphasized is that the dashed line in the figure does 
not represent the laser-driven trajectory averaged over a 
quiver oscillation period but the unperturbed Kepler hy- 
perbola with the identical asymptotic momentum as the 
laser driven trajectory (solid line). Thus, the angular mo- 
mentum of the Kepler hyperbola is identical to that of 
the asymptotic I of the laser-driven electron. Note that 
the number of quiver oscillations along the Kepler orbit 
is not unique thus allowing for path interferences. The 
distance between the hyperbola and the nucleus at the 
pericenter is given by ji^l r m i n = ^^/l + (kl) 2 — ij /k 2 . 
This distance should be equal to the quiver amplitude 
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a. The initial conditions for the 



laser driven trajectory are provided by tunneling ioniza- 
tion with the release of the electron with zero longitu- 
dinal velocity at times U near the maxima of the field 
amplitude F(ti) ~ Fq. Trajectories released at different 
times ti or different maxima of the field reaching the same 
asymptotic branch of the Kepler hyperbola will interfere 
and generate GRT fringes. In order to reach the limit- 
ing case of the dominance in the semiclassical domain of 
a single P/ it is necessary that interference trajectories 
at fixed energy exist that approximately cover the entire 
range of scattering angles (0 < 9 < it) all of which with 
angular momenta close to Iq. Our CTMC-T calculations 
show that close to threshold (k < 0.2) such families of 
trajectories indeed exist. They correspond to Kepler hy- 
perbolas with very similar opening angles of their asymp- 
totes but with the angle of major axis relative to the laser 
polarization, 9, distributed between < 9 < tt. 
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FIG. 3: Partial ionization probability, pi, as a function 
of the angular momentum I for different spectral regions 
indicated in the figure. The parameters of the field are 
lo = 0.05, t = 1005, (a) 7 = 1.34 (F = 0.0377) and (b) 
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7 = 0.67 (F = 0.075). 
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FIG. 4: (a) I distribution of classical trajectories for the 
same parameters as in Fig. 3 (b). (b) Classical trajec- 
tory of laser-driven electron (solid line) and unperturbed 
Kepler hyperbola of same asymptotic E and I (dashed 
lines). Interferences occur for emission at different times 
ti close to different field maxima. 

In conclusion, we have shown that single ionization of 
hydrogen by a moderately strong ultrashort laser pulse, 
in the transition regime from multiphoton to tunneling 
ionization, gives rise to a complex interference pattern 
in the two-dimensional momentum (k z , k p ) plane. While 
at high momenta remnants of ATI rings remain visible, 
at small k near threshold Ramsauer-Townsend diffrac- 
tion oscillations develop at fixed k as a function the an- 
gle between emission direction and laser polarization. A 
simple semiclassical analysis identifies the fringes result- 
ing from interfering paths released at different times but 
reaching the same Kepler asymptote. The present result 
shows that a proper semiclassical description along the 
lines of the "simple man's model" [25| requires a three- 
dimensional description to account for Coulomb scatter- 
ing. Our results feature a striking similarity to recent 
data by Rudenko et al. 0] suggesting the presence of the 



2D interference fringes to be independent of the specific 
atomic core potential. 
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